A New Mean-Field Theory of the Kondo Resonance at Finite Bias 
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We introduce a new slave-boson mean-field treatment of the Kondo effect in a quantum dot 
attached to the leads, when the bias voltage across the leads is finite. The model employs two 
slave boson and two pseudo-fermion operators to express the localized electron. The solution of 
the mean-field equations gives in general two resonance peaks pinned to the chemical potential of 
each lead. The Kondo temperature is shown to scale as min(T^, (Tjf ) 2 /V), where is the Kondo 
temperature at equilibrium, and V is the chemical potential difference of the leads. 



The realization of Kondo resonance in a single molecule 
coupled to metallic leads offers a prime example of the 
many-body correlation effect at work on the nanometer 
scale [1]. The rapid advance in the fabrication technology 
of molecular-scale systems with current-carrying capabil- 
ity will produce growing opportunity to observe Kondo 
effect in a more intricate setup. On the theoretical side, it 
is important to develop tools suitable to describe Kondo 
phenomena in situations of increasing complexity. 

The Kondo effect in a quantum dot coupled to metallic 
reservoir has been thoroughly studied, and many theoret- 
ical tools are available [2] . The formation of a resonance 
state at low temperature (Abrikosov-Suhl resonance) at 
the Fermi level of the host metal can be understood as a 
consequence of hybridization between a localized electron 
and the host electrons. The slave-boson theory captures 
many of the essential features of the Kondo resonance in 
the so-called Fermi-liquid regime [3] . For a quantum dot 
coupled to leads with different chemical potentials, i.e. 
Hl Mi? where fir, and /j,r are chemical potentials of 
the left and the right leads, however, there is a discrep- 
ancy in the outcome of the slave-boson mean-field the- 
ory(SBMFT) (predicting a single Kondo resonance) [4] 
and that of equation-of-motion(EOM) or non-crossing- 
approximation(NCA) theories (predicting a splitting of 
Kondo resonances) [2]. 

It is important to resolve this discrepancy fully, in view 
of the simplicity (hence its utility) of SBMFT with which 
to capture essential features of the Kondo physics. It will 
be desirable if the outcome of the mean-field theory can 
be made to agree with that of the more advanced tech- 
niques such as NCA or EOM methods. In this paper, we 
assume that the predictions of the NCA/EOM is indeed 
correct [5] , and offer a new version of SBMFT that gives 
rise to split Kondo peaks. We conclude, on the basis of 
the new mean-field theory, that (I) the spectral function 
of the quantum dot has resonances pinned at the chem- 
ical potentials of each lead, (2) the Kondo temperature 
scales as min(T^, {T^) 2 /V) where T£ is the Kondo tem- 
perature in the case of un-biased leads coupled to the dot, 
and V is the chemical potential difference of the leads. 
In the following we present the mathematical formula- 
tion of the theory and discuss physical consequences of 
the solution that emerges. 



A quantum dot coupled to two leads via tunnelling 
may be modelled by the Hamiltonian 

H = H L +H R + H d + H T (1) 
where the lead, the dot, and the tunnelling parts are 
h l(R) = E ( £fc ~ VL(R))ct [R)klT c L{R)k!y 

k£L(R),a 



H T = V L 



E 

k£L,a 



d+d a + Un d (n d - I) 
V R 



r + d 



E 

k£R,a 



c+ ka d a + h.c. (2) 



Tunnelling amplitude to each lead is given by Vl and Vr 
while the density of states of leads are assumed to be 
constant, pi, and p R respectively. We consider the case 
where the bare energy level of the dot — e d (< 0) and the 
energy of double occupancy U — e d (> 0) are both large 
(Kondo regime). 

The d-level operator d£ in the U — *• oo limit is really 
rf+ = |er)(e|, taking an empty state (|e)) into a singly 
occupied state of spin a ( because doubly-occupied 
state is energetically forbidden. Commutation relations 
for d a thus differ from the canonical one: 

K,d+} = l-tf± CT <L CT , {d a ,d+_ a } =d+_ a d a . (3) 

In the usual slave-boson approach, the localized electron 
operator d a is replaced by b + f a with b and f a being 
the standard boson and fermion operators. Commuta- 
tion relations, Eq. (3), are recovered provided we require 
b + b+^2 a fa = 1- The d a operators in the Hamiltonian 
are replaced by an equivalent expression, d a — > b + f a , 
while the constraint is included in the form of a Lagrange 
multiplier, 



H 



H' = H + \(b+b + J2faU-l) 



(4) 



The free parameters A and b are self-consistently deter- 
mined by the condition 



d(H') 



= = 



d(H') 
db 



(5) 



The stationary state solution describes a narrow reso- 
nance state formed at the chemical potential of the lead 
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and gives a quantitative estimate of the width of the reso- 
nance that is also associated with the Kondo temperature 

1 K- 

The same approach was adapted by Aguado and Lan- 
greth (AL) to the non-equilibrium case, ^ Hr [4]. 
One finds in this case that a single Kondo resonance state 
is formed at the energy level between the chemical po- 
tentials of each lead. In physical terms, the chemical 
potentials of the leads are "averaged over" , before an in- 
teraction with the localized electron is considered. 

The root of the Abrikosov-Suhl resonance is the high- 
order virtual hopping process that in effect exchanges 
spins between the lead electrons and the localized elec- 
tron. When more than one lead is attached to a dot, 
it appears that a correct picture is to view such hopping 
processes as occurring independently with respect to each 
lead, as far as one can ignore a direct interaction between 
the leads. Naturally, then, the resonance state will have 
to be formed for each lead, at the chemical potential of 
the corresponding lead. Motivated by this picture, let us 
write the ci- level electron operator as 



(6) 



Two sets of slave-boson (61,62) and pseudo-fermion 
(/i<7,/2ct) operators have been introduced, with relative 
weights a and (3 to be determined shortly [6]. The com- 
mutation algebra, Eq. (3), are recovered provided we 
require the following: 

(a) 6+6 1+ ^/+/ la =p 

<7 

(6) b+b 2 + J2.flf2« = q 



(c) dtd a = a 2 f+J la + p 2 f+J 2a 

(d) a 2 p + f3 2 q=l. 



(7) 



Two positive constants p and q have been introduced in 
the above relations. The empty state \e) takes on the 
representation, 



|e) = (m6+6+ + ^(6+) 2 )|0), 



(8) 



where |0) is the vacuum state for the slave boson and the 
pseudo-fermion operators, and for the singly-occupied 
state, |cr) = d+\e). This assignment is very different 
from |e) = 6 + |0) and \a) = f^\0) of the conventional 
slave-boson theory. Note also that the expression for the 
empty state is asymmetric in the interchange of 61 and 
62 [7]. Later we will see that this asymmetry is a conse- 
quence of the finite bias across the dot. 



In order to ensure that the new Hamiltonian expressed 
in terms of slave bosons and pseudo-fermions has the 
same matrix elements as the original H, we further re- 
quire that (er|<i+|e) = 1 in the new representation. In 
meeting this requirement, and (a)-(b) of Eq. (7), we 
deduce several relations: (1) p = 1 + k, q = 1 — k, 



(2) 



= 1 - k, n 2 = k, (3) a 2 = /3 2 = 1/2. Now 



the d-level operator in Eq. (6) is uniquely defined as 
d a = {b~l fia + 6j/ 2(T )/v / 2. The value of k will be deter- 
mined later. 

Hamiltonian in the slave-boson representation becomes 
Ht = \ E c+ LkMh° + btf2*) + h.c. 

V keL,a 



+ \ E c RkMha + b+ha) + h.c, (9) 



V2 



keR,a 



and 



(10) 



The complete Hamiltonian together with the occupation 
constraints given in Eq. (7)(a)-(b) becomes 

H = H L + H R + H d + H T 

+ Ai ^6+6!+^/+/^-!-^ 



+ A 2 



6+62 



(11) 



In the following we absorb \[2 in the definition of the 
tunnelling amplitudes, Vl and Vr. 

We have accomplished an exact mapping of the origi- 
nal Hamiltonian, Eq. (2), to an equivalent one, Eq. (11), 
using the slave-boson theory valid in U — > 00 limit. Al- 
though the one-slave-boson theory of AL also starts from 
an exact mapping of the original Hamiltonian [4] , the out- 
come of the mean-field analysis may differ significantly 
depending on the type of mapping chosen. 

One can carry out the stationary-state analysis of the 
ground state of the above Hamiltonian in the manner of 
AL, solving for the four parameters, 61, 62 and Ai and 
A2, which requires 



and 



b i(2) + (E /l(2)a/l(2)a) = 1 ± K, 



2Ai(2)6 1 ( 2 ) + (H T i( 2 )) = 0. 



(12) 



(13) 



The tunnelling part of the Hamiltonian factorizcs, Ht = Hti + Ht2, with the 1 and 2 parts containing the 1- or 
2-dependent portion of the d-level operator. We replaced the 6-operators by their mean-field values 61 and 62 that 
are assumed to be real. Keldysh techniques are employed in evaluating the averages [2,4]. At zero temperature, the 
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mean-field conditions of Eqs. (12) and (13) turn into 



2 v 2 ' T L +l 



■ arctan 



■ arctan 
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A 2 - Pl 
6 2 £ 2 
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R 



Ai - pr 

A2 — PR 
b^2 



(14) 



and 



e d = T L In 



e d = T L In 



D 2 



(Ai - ml) 2 + (& 2 S 2 ) 2 
D 2 

(A 2 - /i L ) 2 + (fo 2 £ 2 ) 2 



r fl in 



+ r fl in 



D 2 



(Ai - p R ) 2 + (&?S 2 ) 2 
£> 2 

(A 2 - M K) 2 + (fo|S 2 ) 2 



(15) 



We have introduced Tl = V 2 Pl,^r = VrPr, the width of the conduction band D, the self-energy E 2 = k{Yl + r#), 
and re-defined Ai( 2 ) — e<j + ^1(2)^1 as Ai(2) in the above expressions. £ = Si + z£ 2 gives the complex self-energy of 
the (i-electron due to hybridization with the leads. We consider the case e d 3> Ai( 2 ), |S|, as would be appropriate for 
the Kondo regime. 



Equations (14)-(15) are the central results of this pa- 
per. In the following we solve them for the case of sym- 
metric coupling, Tl — Tr, and in the limit of infinite 
barrier to one of the leads, Tr — ► 0. 

The bias voltage pr — pl = V is taken to be pos- 
itive (right lead at higher chemical potential than the 
left lead). We look for solutions in the Kondo limit 
b\ w w b 2 , and Ai and A 2 lying very close to either 
PL or pr. The case Ai w A 2 has a solution only for a 
small bias, V <C b 2 Y>2 an d is n °t applicable at large V. 
Hence we require that Ai w /Ul, and A 2 w /Ur. The other 
possibility, Ai w and A 2 w /Uj,, is obtained when the 
bias is reversed, V < 0. 

The self-consistency equations are now reduced to 
(taking T L = T R ) 



b\^bl/[l + , 2 {bl-n') 2 Yl 2 
b 2 /[l+n 2 (b 2 + n') 2 } 1/2 , 



where «/ = n — arctan[V r /6§S 2 ], and 

D 2 



b 2 Z 



0^2 



-ircd/S 2 



maa^F, &gS 2 ) 
while Ai — pl and A 2 — satisfy 

Al(2) - PL(R) ~ ^1(2)^2 x fe 2 (2) - k' 



(17) 



(18) 



Now we come to the determination of n. It is cho- 
sen to minimize the energy (Ht + H d ) which, after self- 
consistent relations Eqs. (12) and (13) are taken into 
account, become — ^e d (b 2 + b 2 ) + const. Other terms 
in the expectation value of the Hamiltonian are insensi- 
tive to the choice of k. Optimal k is therefore one that 
maximizes b\ + b\. Given the expressions for b\ and b\ 
obtained in Eq. (16), it immediately follows that k' must 
be zero, otk = arctan[y/& 2 £ 2 ]. Since the empty state 



|e) is defined as — nb^b^ + y/n/2(b^ ) 2 )|0), we see 
that the definition of the state itself depends on the ap- 
plied bias! Put another way our expression for the empty 
state contains a variational parameter k, which is chosen 
to maximize the energy gain from the hybridization. 

The mean-field solution demonstrates the existence of 
two Kondo peaks since the <i-electron spectral function 
is given by the sum 



h 2 



(19) 



with the spectral functions for the pseudo-fermions fi 
(Ai(u>)) and / 2 (A 2 (u)) each producing a resonance at 
the chemical potential of the lead [8] . 

The magnitude of 6qS 2 is the so-called Kondo temper- 
ature, Tk- One can see from Eq. (17) that Tk changes as 
the bias voltage is increased through the threshold value, 



(16) V = L> e -^/ 2S 2 



T K in the manner 



•k 



-K 



(T K f/V 



(V < T K ) 
(V > T* K ). 



(20) 



Coleman et al. [9] made a similar prediction for the 
crossover behavior of the Kondo scale as the bias voltage 
is increased through T K , based on the perturbative cal- 
culation of the magnetic susceptibility. Our mean-field 
calculation demonstrates the crossover behavior explic- 
itly. The crossover scale T K is the Kondo temperature 
at equilibrium. 

When tunnelling to one of the barriers is forbidden, 



0, we have as the solution 



0. 



H^2 



T. 



K> 



and Ai = A 2 w pl- Therefore, the two spectral functions 
A\(oj) and A 2 {uo) merge to produce a single resonance 
peak at pl, as one would expect in the case of coupling 
to a single reservoir. Further calculation (see next para- 
graph) shows that the split resonance positions are close 
to the chemical potentials of each reservoir as long as the 
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tunnelling strengths, given by Tl and Tr, are reasonably 
close to each other. As Tr/Yl becomes too small, even- 
tually the resonance level pinned to /x# begins to float, 
and merge with 

When V > T K , some understanding can be achieved 
for the slightly asymmetric case, Yl+Yr = Y, Yl — Yr = 
Ye with < e <C 1. Here we obtain that the Kondo 
temperatures scale differently for the two resonances, 
with the Tk for the stronger resonance near /xl vary- 
ing as T^(T^/V')( 1 - e )/( 1+e \ and for the weaker reso- 
nance to the right reservoir, as T£(T£/U) (1+£)/(1 ~ c) . 
In the highly asymmetric case, i.e. Yt, = Y(l — e) 
and Yr = Te, we get the Kondo temperature scaling 
T£(T£/I0 e/(1 ~ £) for both resonances. Predicting such 
dependence of power-law exponents on the tunnelling ge- 
ometry would be very difficult based on other calculation 
schemes, and illustrates the utility of SBMFT. It is of 
course an interesting issue to see if experimentally such 
power-law behavior indeed shows up. 

In conclusion, we have introduced a new version of 
SBMFT for the Kondo effect in a quantum dot that cor- 
rectly reproduce the equilibrium Kondo effect, and inter- 
polates smoothly to the non-equilibrium situation with a 
bias V or to the case of asymmetric coupling to the reser- 
voirs. The introduction of two sets of slave-bosons and 
pseudo-fermions may at first sight seem to have doubled 
the number of orbitals in the dot, but on closer inspec- 
tion it is clear that the dimension of the Hilbert space 
for the dot is strictly three, as ought to be the case for 
U — > oo. Rather, we believe this doubling of the number 
of operators is tied to the presence of multiple reservoirs 
with which the dot's electron can interact. Meanwhile, 
the technical aspect of the present theory is completely 
straightforward and may be easily generalized to under- 
stand Kondo physics in other, more complex situations 
with multiple dots [4], a phonon-coupled dot [10], etc. 
Also of great interest is the role of the inelastic processes, 
not considered in the present mean-field theory, in reduc- 
ing the Kondo scale Tk [11]. Future work will address 
these issues. 
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